Intermittent distribution of tracers advected by a compressible random flow 
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Multifractal properties of a tracer density passively advected by a compressible random velocity 
field are characterized. A relationship is established between the statistical properties of mass on 
the dynamical fractal attractor towards which the trajectories converge and large deviations of the 
stretching rates of the flow. In the framework of the compressible Kraichnan model, this result is 
illustrated by analytical calculations and confirmed by numerical simulations. 



We are interested in the passive transport of a scalar 
density by smooth-in-space compressible random flows 
in a d-dimensional bounded domain A. During trans- 
port, the density develops strong inhomogeneities. Our 
goal is to describe quantitatively their fine structure aris- 
ing at asymptotically large times. Compressible flows are 
physically relevant not only at large Mach numbers. For 
instance, the dynamics of a suspension of finite-size (iner- 
tial) particles in an incompressible flow may be approxi- 
mated in the limit of short Stokes times bythat of simple 
tracers in an effective compressible flow 0, 0] ■ Another 
application is to the advection of particles floating on the 
surface of an incompressible fluid 

A smooth passive density field p(t, x) evolves in the 
velocity field v(t, x) according to the continuity equation 



d t p + V • (pv) = 



(1) 



which preserves the total mass that we shall assume for 
convenience normalized to one. If the initial density at 
time io is uniform and if the flow is compressible and 
sufficiently mixing, then the solution p(t) of Q will ap- 
proach, when to tends to — oo, a singular limit p*(t) which 
is a measure with support on the dynamical attractor to- 
wards which the Lagrangian trajectories converge. For 
random velocities, the measure p*(t) and the dynamical 
attractor depend on the velocity realization. We shall 
consider stationary velocity ensembles where the statis- 
tics of p*(t) does not depend on t. For convenience, we 
shall restrict our study to = p*(0). One expects the 
measure p, to have roughly a local product structure 
with a smooth density along the unstable manifolds of 
the flow and a fractal-like structure in the transverse di- 
rections (such measures are called SRB Q). 
The Lagrangian average defined by 



(F) = / F(x\v)p*(x\v)dx, 

JA 

where the overline denotes the expectation with respect 
to the velocity ensemble, samples points in A according 
to the density p» of the asymptotic tracer distribution on 
the random attractor. It determines an invariant mea- 
sure /i* of the random dynamical system defined on the 
product of the physical space A and the space of the ve- 
locity realizations. 



We are interested in the small-scale statistics of the 
mass distribution associated to the measure p*. Denote 
by B r (x) the ball of radius r around the point x and 
introduce the quantities 



m r (x) = / p*(y)dy , 

JB r {x) 



In m r (x) 

hr{x) = : — 

mr 



The local dimension at x, defined if the limit exists as 
h(x) = lim r ^o h r (x), characterizes the small-scale dis- 
tribution of mass associated to the limiting density p*. 
A more global assessment is provided by the Hentschel- 
Procaccia spectrum for dimensions (a, ||| of p* : 



HP{n\p*) = lim 



ln/ A m r (x) n 1 p*(x)dx 
(n — 1) lnr 



(2) 



for real n. In particular, HP(0 \ p*) is the fractal dimen- 
sion of the support of p x , HP(1 \ p*) is known as the in- 
formation (or capacity) dimension and HP (2 | p„) as the 
correlation dimension. For random flows, one may de- 
fine the quenched version of t he Hcntschcl-Procaccia 
spectrum HP qu (n) = HP(n | p*) and the annealed ver- 
sion HP an (n) where the velocity average is taken under 
the logarithm in the numerator on the right-hand side 
of ©. One has (n - l)HP qu (n) > (n - l)#P an (n) 
and _ff_P qu (f) = HP an (l) (if exists). If the moments of 
m r (x | v) with respect to the invariant measure /x* ex- 
hibit the small r scaling 



(3) 



then HP an (n + 1) = £n/n. There are various aspects of 
behavior ©. First, a non-linear dependence of the scal- 
ing exponents £ n implies intermittency in the mass dis- 
tribution and, in particular, the presence of non-Gaussian 
tails in the probability density function (PDF) of m r . 
Second, this behavior suggests that the PDFs of h r take, 
for small radii, the large deviations form e( lnr ) s (' 1 ) with 
the rate function S(h) and the scaling exponents re- 
lated by the Legendre transform: S(h) — max,, [£« — nh]. 
Third, the (annealed) fractal dimensions f(h) of the ran- 
dom sets on which the local dimension h(x) is equal to h 
(the multifractal spectrum), is expected to be equal 
to h - S(h). 
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The aim of this letter is to relate quantitatively the 
dimensional spectrum HP an (n) to the large deviations of 
the stretching rates of the flow (sometimes also called 
finite-time Lyapunov exponents). We note here that in 
a similar relation was discussed for n = 1 and, in a 
qualitative way, for large n. 

Let X(t,x) denote the Lagrangian trajectory passing 
at time zero through point x. The matrix W(t,x) = 
(djX l (t, a;)) describes the flow linearized around the tra- 
jectory X. The eigenvalues of the positive matrix W T W 
may be written, arranged in non-increasing order, as 
e 2t<Tl , . . . , e 2tt7d with G{ = (Ti(r, x\v) called the stretch- 
ing rates of the flow. The limits 

Xi = lim <7i(i) 

t — >oo 

exist almost surely with respect to the invariant measure 
//* and are constant if the latter is ergodic. They define 
the Lyapunov exponents of the flow. This is the content 
of the Multiplicative Ergodic Theorem 0. As argued 
in p|, if the exponents A^ are all different then the PDF 
of the stretching rates takes the large deviation form 



-t H(a 1 ,..., a d ) 



((71 - (72 ) 



(4) 



for large t with convex rate function H attaining its min- 
imum, equal to zero, at cr^ = A^. 

If the top Lyapunov exponent Ai is negative, then the 
attractor measure p» is expected to be atomic with triv- 
ial mass statistics. For Ai > 0, Ledrappier and Young 
showed under rather general assumptions that the lo- 
cal dimension h is almost surely equal to the Lyapunov 
dimension = j + 5 where < S < 1 is such that 
Ai + . . . + Xj + <5Aj+i = 0. The above statement implies 
that 



«P qu (l) = aP an (l) = d£„/dn|. 



(5) 



The Lyapunov dimension was introduced by Kaplan 
and Yorke 10] and may be heuristically interpreted as the 
dimension of objects keeping a constant volume during 
the time evolution. Eq. 10 gives only partial information 
about the mass scaling as compared to the full set of 
exponents £ n . 

To relate the mass distribution in small balls to the 
fluctuations of the stretching rates, we focus on the two- 
dimensional case. The expected smoothness of in the 

time time t 




FIG. 1: Two-dimensional sketch of the backward-in-time evo- 
lution between times t and of a small parallelogram of length 
ri in the direction of the unstable manifold and T2 in the di- 
rection perpendicular to it. 



unstable direction leads us to consider small parallelo- 
grams with one side of size r\ parallel to the unstable 
manifold and with extension r 2 in the direction perpen- 
dicular to it. We expect that the moments of the mass 



»ri,r 2 



in such parallelograms obey 



(G) 



for small ri, r 2 as long as £„ — n > 0. When < n, we 
expect © to be replaced for r 2 > T\ by 

(<U)~ri"r 2 °. (7) 

This may be viewed as analogous to the stretching in one 
direction of a fractal set of dimension D, while contract- 
ing it in the other one. In the expanding direction, the 
set will behave (down to a scale depending on the stretch- 
ing) as if projected on a line, so it will have dimension 
min(l, D). In the other direction it will have the comple- 
mentary dimension D — min(l, D) — max(Z) — 1, 0). 

Let us consider such a parallelogram at t > and let us 
look at its pre-image at time zero. While the direction 
parallel to the unstable manifold is exponentially con- 
tracted backward-in-time with a rate given by the largest 
stretching rate a\ , the other direction typically expands 
with an exponential rate <J2- Hence the time-zero pre- 
image of the original parallelogram is (approximately) 
another parallelogram as sketched in Fig. Conserva- 
tion of mass and stationarity of the statistics lead to the 
relation 

If there is a sufficiently rapid loss of memory in the La- 
grangian dynamics (i.e. if the invariant measure is suf- 
ficiently mixing) then the expectation on the right-hand 
side should factorize for large t (such a factorization holds 
for all t in the Kraichnan model discussed below). In such 
situation, using the large deviation form of the PDF 
of the stretching rates, we infer that 

(< ie -w, r3 e-^)e" m(CT1 '-' CT2) a<nd* 2 . 

Consistency of the above relation with the scaling © 
requires that 

1- j e~ t[mTl+( ^-" )CT2+ff(,Tl ' <J2)1 deride • 

J <T1 >0"2 

Since t is assumed large, a saddle-point argument implies 
then the following relation between the scaling exponents 
and the rate function H of the stretching rates: 

min \no\ + (£„ - n)a 2 + H(ai, a 2 )} = 0. (8) 

ax >o"2 

Alternative formulations are 

£,n = n — max — [nai + H(ai, a 2 )] — min [hn + S(h)] 

o"l>C2 (72 ' l >0 

O2<0 
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for S(h) = min (T> o a^ 1 H((h — l)rr, —a). These formulae 
are valid for £ n > n. Similarly, from we may get the 
formula valid for £„ < n 

£ n = - max —H(ai,a 2 ). (9) 



max 

(71 ><T2 
<T1<0 



01 



It is easily checked from |jH)l that (d£ n /dn)|„ = o = 1 — 
A1/A2 which coincides, in our settings, with formula JSJ. 
Also from (JSJ one may deduce that 



lim ~{\R(t)\- 



0, 



for R(t) — W(t, x)Ro, meaning that the generalized Lya- 
punov exponent of order — £1 vanishes. This is a known 
result for stochastic flows of diffeomorphisms, see [Tl) . 
Relation |JSJ can be extended to dimensions higher than 
two using similar arguments. 

To illustrate the small-scale properties of the mass dis- 
tribution, we focus in the sequel on transport by a ran- 
dom compressible velocity field chosen in the framework 
of the Kraichnan model [T^] • For the domain A we take a 
periodic box. The velocity v is taken centered Gaussian, 
with covariance 



v l (x + £,t + T )v3(x,t) = 2(IV y -tf°{e)) 5(t), 
where, for small separations £, the function d lJ satisfies 

= Hi + 1 - 2p) 8 l H 2 + 2{pd - 1) tl j + o(£ 2 )] 

assuring local isotropy. The parameter p, called the com- 
pressibility degree, is chosen in the interval [0, 1]. Its ex- 
treme values correspond, respectively, to incompressible 
and to potential velocity fields. For such velocity ensem- 
bles, the distribution of the stretching rates is explicitly 
known 0, ^| . The corresponding rate function takes the 
simple form 



H = 



i-[]r>-A,) 2 + C 2 ( 

i— 1 



d 

E 

i=l 



d = 4Di (d+p(d-2)), C 2 ee (l-pd)/(p(d-l) (d+2)). 
The Lyapunov exponents are Xj = Di[d(d — 2j + 1) — 
2p(d+(d— 2)j)]. In the two-dimensional case, the largest 
Lyapunov exponent is negative when p > 1/2. For p < 
1/2, the exponents for the mass distribution obtained 
from © and © read 



2rt+ A /(l+2p) 2 -8pn 



1 



if 
if 



n < n r 



n > n c 



(10) 



where the critical moment n cr and the saturation expo- 
nent are given by 



1 



2 V 2p 
00 = 2?T- cr 1 



if 0<p<l/6, 



e- = ei = f^ if i/6<p<^ 



The two different behaviors are illustrated in Fig. 2. In 
both cases, the events contributing to the saturation of 
the exponents are those for which a mass of order unity 
is concentrated inside the small ball. 
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FIG. 2: (a) Ellipses in the (n, f„) plane corresponding to 
the unconstrained minimum in Dotted parts refer to 

situations when either the minimum is reached for o\ < 02 or 
< n. Two values of p are represented to illustrate both the 
case p < 1/6 and p > 1/6. (b) Location where the minimum 
is reached in the (<ti, 02) plane. 

Numerical simulations confirm the values of the scal- 
ing exponents obtained explicitly for the compressible 
Kraichnan model. To distinguish the two cases, two 
different values of the compressibility degree are inves- 
tigated (p = 1/10 < 1/6 and p = 3/10 > 1/6). The 
velocity field v is generated by the superposition of nine 
independent Gaussian modes and the density is approxi- 
mated by considering a large number of Lagrangian trac- 
ers. The exponents obtained numerically after averaging 
over 10 5 turnover times and for N = 10 5 tracers are 
shown in Fig. |3J Although statistical convergence of the 
average is quite slow, these expensive simulations are in a 
rather good agreement with the theory, in particular with 
the saturation of the exponents after the critical order. 





(a) p = 1/10 



(b) p = 3/10 



FIG. 3: Scaling exponents £ n for the mass distribution as- 
sociated to the advection of a density field by a Kraichnan 
velocity field. For two different degrees of compressibility, the 
exponents obtained match those predicted by theory which 
are represented as dashed lines. 

The determination of the exponents may be improved 
for positive integer orders n by considering only n+1 par- 
ticles. This method allows to perform very long time av- 
erages (here of the order of 10 8 turnover times) required 
for good convergence of the statistics at small scales. As 
shown in Fig.^J for p = 3/10, the moments clearly scale 
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over several decades. For the lower value of compress- 
ibility p = 1/10, convergence is slower because of the 
smaller probability for the tracked particles to come close 
together. 



(a) (1=1/10 (b)p = 3/10 




FIG. 4: For n + 1 particles (n = 1, 2 and 3) the distribution 
of the maximum of the distances of the reference particle to 
the other n particles is represented. It can be shown that this 
PDF is proportional to d(m")/dlnr. The insets show the 
difference with small-r asymptotes with slopes given by the 
theoretical exponents. 

Finally, let us mention that, for integer orders, the 
exponent £ n for the Kraichnan model can be linked to the 
homogeneity degree of the stationary single-time density 
correlation function 



Fn+i = p(0)p(xi) ■ ■ ■ p{x n ) . 
The n th order moment of m r may be written as 

(to™) = / .F n +i(0, xi, . . . ,x n ) dxi . . .dx n . (11) 

J Z3 r (0)x---xB r (0) 

In the stationary regime of the Kraichnan model T n+ i is 
a zero mode of the operator 

M Li = E n ■.'>■, [( D ° sij - dij ( x » - ^)) • ] 

0<k,i<n 

We need not write these zero modes explicitly. The ho- 
mogeneity degree of the isotropic solution of lowest de- 



gree can be found simply by requiring its positivity and 
imposing on it certain continuity and integrability con- 
ditions. The branch n < n cr in (|10[1 is obtained by re- 
quiring the solution to be continuous at x\ = . . . = x n . 
On this subspace, M^ +1 is degenerate, owing to the fact 
that collinear points remain collinear when transported 
by the (linearized) flow. If in addition we restrict M^ +1 
to the rotationally invariant sector, it becomes an or- 
dinary second-order homogenous differential operator. It 
has two scaling solutions. The one with exponent £ n — nd 
gives, through l|ll|) . the non-saturated (n < n cr ) branch 
of 1 11 (J| l. Such a solution breaks down when the corre- 
sponding zero mode ceases to be integrable at small Xi 
which may occur if the restriction of the zero mode to the 
collinear sector is not integrable near 0, i.e. — nd < — 1. 
This gives n cr and £00 in the case p < 1/6. Another sce- 
nario is if the zero mode has a non-integrable singularity 
around the collinear geometry. This gives n cr and ^ in 
the case p > 1/6. 

An important open question not touched upon by this 
paper concerns the case of rough-in-space velocity helds 
appearing in the limit of very high Reynolds numbers. 
That problem cannot be formulated in terms of stretch- 
ing rates but the relationship with density correlations 
and zero modes still holds. 
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When this work was essentially finished, we learned 
from A. Fouxon that he has also derived the relationship 
© between the scaling exponents of mass and the large 
deviations of the stretching rates. 
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